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pLn ' Abstract. The dual origin of population II Li, in both big bang nu- 

Q . cleosynthesis and galactic cosmic-ray nucleosynthesis is discussed. It is 

psj ' argued that with additional ^Li data, stringent limits on the degree of 

^Li depletion can be obtained. ^Li depletion is also constrained by the 

concordance of big bang predictions with observational determinations 

of light element abundances. Stringent limits can also be obtained for a 
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1. Introduction 



The key link between galactic cosmic-ray and big bang nucleosynthesis (GCRN 

and BBN) is the production of ^Li. Although the big bang produces each of 

O \ the element isotopes ^Li, ^Li, ^Be, ^'^B, and ^^B, their abundances relative to 

GCRN production are nearly negligible, with the exception of ^Li for which 
BBN is the main source for the abundances determined from the observation of 
population II halo stars. In contrast, the other LiBeB elements are all formed 
predominantly in GCRN. Therefore, not only is ^Li a common link relating the 
^ ' two process, but the interpretation of the abundance of ^Li in halo stars may 

require an understanding of both mechanisms. 

One of the goals of this paper is to differentiate the sources of ^Li and in 
particular the primordial abundance of ^Li. Among the criteria that will be used 
is the concordance of the primordial ^Li abundance with the other light element 
data (^He and D), as well as the consistency with the galactic evolution of ®Li 
and BeB. BBN produces significant amounts of D, ^He, ^He, and ^Li (see e.g. 
Olive 1999). Furthermore, in the standard BBN model there is really only one 
free (undetermined) parameter, the baryon-to-photon ratio, r]. Therefore, we 
can impose as a constraint the concordance between the predicted abundances 
of the light elements with the observationally determined abundances (with the 
exception of '^He which is highly dependent on the uncertainties of stellar and 
galactic evolution). Li, however, is also produced in spallation and fusion pro- 
cesses associated with cosmic rays. By modeling the evolution of the observed 
abundances of Be and B, we can determine the associated ^Li production (in a 
given model). As we will see, ^Li has the potential to play a key role in this type 
of investigation, since its GCRN production is very similar to that of ^Li and is 
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much less model dependent. However, the current paucity of data (though we 
should celebrate the fact there is some data) makes it difficult at this time to 
draw hard conclusions. 

In what follows, we will very briefly review the key essentials of BBN rele- 
vant to our discussion here. We will then discuss the relevant LiBeB data and 
the dual origin for ^Li. The main focus of the paper will be devoted to the 
question of ^Li depletion and the limits that can be placed on depletion from 
BBN and GCRN. 



2. Big Bang Nucleosynthesis 

As noted above, the standard BBN model is in fact a one-parameter model. 
The key uncertainty in the prediction of the light element abundances is the the 
baryon-to-photon ratio, rj. The other parameters normally associated with BBN 
are now to a large extent fixed. Experiments at LEP have fixed the number of 
light neutrino fiavors, Ni, = 3, and the neutron mean- life is now well measured 
to be Tn = 886.7 ± 1.9 s (see e.g. the Review of Particle Properties, 1998). The 
predicted abundances of the light elements are shown in Figure 1, which concen- 
trates on the range in r/io between 1 and 10 (ryio = W^^rj). The curves for the 
^He mass fraction, Y, bracket the computed range based primarily on the uncer- 
tainty of the neutron mean-life. Uncertainties in the produced ^Li abundances 
have been adopted from the results in Hata et al. (1996). Uncertainties in D and 
^He production are small on the scale of this figure. The boxes correspond to 
the observed abundances with 2a statistical uncertainties and the dashed boxes 
include systematic uncertainties. 

As one can see, the predicted abundances of the light elements are roughly 
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Yp ^ 0.23 - 0.25 

D/H Ki 3.4 - 20 X 10~^ 

*He/H ^ 1 - 2 X 10"^ 

'Li/H « 1 - 2 X 10"^° (1) 



To test the concordance of the BBN predictions with the observations it has 
become useful to use a maximum likelihood method as described in detail in 
Fields et al. (1996). Shown in Figure 2a, are individual likelihood distributions 
for ^He and ^Li. To obtain these distributions, the observed values of Yp = 
0.238 ± 0.002 ± 0.005 (Olive, Skillman, & Steigman 1997, Fields & Olive 1998) 
and ^Li/H = (1.6 ± 0.1) x 10"^° (Molaro, Primas, & Bonifacio 1995, Bonifacio 
& Molaro 1997) were adopted. ^Li exhibits a two-peaked distribution due to the 
local minimum in the predicted abundance versus rj. These two elements are 
shown because of their relative lack of model dependent galactic evolutionary 
effects, in contrast to D and ^He. In Figure 2b, the total likelihood function 
is shown (the product of the two distributions). As one can see there is broad 
agreement (based on ^He and ^Li) in the range 1.55 < 7710 < 4.45. 

In Figure 3a, the total likelihood distribution including a high value of 
deuterium D/H = (2.0 ±0.5) x 10~^ is shown (Songaila et al. 1994, Carswell et 
al. 1994, Webb et al. 1997, Tytler et al. 1999). Now, the concordance is limited 
to a narrower range in rj, 1.5 < rjiQ < 3.4 as can be seen from the total likelihood 
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Figure 1. The light element abundances from big bang nucleosyn- 
thesis as a function of rjiQ = W^^t]. 



distribution shown in Figure 3b. Similarly, in Figure 4, the distribution is shown 
for D/H = 3.4 ±0.3 x lO'^ (Buries & Tytler, 1998a,b). Notice the relative scale 
in the two distributions (3b and 4b). The low value of D/H is compatible with 
He and ^Li at the 2 a level in the range 4.2 < 7710 < 5.6. 

In contrast to the isotopes discussed above, big bang nucleosynthesis also 
produces some ^Li, Be, and B. In the favored range of 1.5 - 4.5 for r^io, the big 
bang produces roughly (Thomas et al. 1993), 



^Li/H 
^Be/H 
ioB/H 
"B/H 



2 - 9 X 10^^^ 
0.04 - 2 X 10"^"^ 
0.5-3 X 10"^^ 
0.02 - 1 X 10"^^ 



(2) 



Because these abundances are far below the observed abundances found in 
Pop II halo stars, (^Li/H ^ few x 10"^^ ^Be/H ~ 1 - 10 x 10-^^ and B/H 
~ 1 — 10 X 10"^'^), it is generally recognized that these isotopes are not of primor- 



dial origin, but rather have been produced in the Galaxy, through cosmic-ray 
nucleosynthesis . 





Figure 2. (a) - Likelihood distribution for each of ^He and ^Li, 

shown as a function of rj. The one-peak structure of the ^He curve 
corresponds to the monotonic increase of Yp with rj, while the two 
peaks for ^Li arise from the minimum in the ^Li abundance prediction; 
(b) - Combined likelihood for simultaneously fitting ^He and ^Li, as a 
function of i]. 





Figure 3. (a) - As in Figure 2a, with the addition of the likelihood 
distribution for D/H assuming high D/H; (b) - The total likelihood 
distribution with high D/H. 



3. LiBeB Data 

The LiBeB data in population II halo stars is been discussed extensively in these 
proceedings and therefore, we can be brief here and concentrate on the data as 
it concerns nucleosynthesis and the question of the primordial abundance of ^Li. 
The Li data was reviewed by Molaro, and based on the work of Molaro, 
Primas & Bonifacio (1995) and Bonifacio & Molaro (1997) as well as more recent 
results, he argued for a uniform abundance of ^Li with a value 



^Li/H~ (1.6 ±0.1) 



10 



-10 



(3) 



absent of any trends with respect to either temperature or metallicity ([Fe/H]). 

The ^Li abundance (reviewed here by Hobbs) has been improving over time. 

There are now two stars at low metallicity ([Fe/H] ~ —2.3) with measured ^Li/H 
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Figure 4. (a) - As in Figure 2a, with the addition of the hkehhood 
distribution for D/H assuming low D/H; (b) - The total likelihood with 
low D/H. The dashed curve is the distribution from Figure 3b. 



(Smith, Lambert & Nissen 1993,1998, Hobbs & Thorburn 1994,1997, Cayrel et 
al. 1999): 



HD84937 : 
BD26°3578 : 



5Li/Litotai = 0.054 ±0.011 
^Li/Litotai = 0.05 ± 0.03 



(4) 



corresponding to ^Li/H ~ 6 — 10 x 10^^^. These values, together with the solar 
abundance (^Li ~ 1.5 x lO^^'^) lead to strong constraints on the evolution of 
^Li and more importantly on the nucleosynthesis of the LiBeB elements and the 
depletion of ^Li as will be discussed below. 

The BeB data was discussed by Duncan and Garci-Lopez. The pop II BeB 
data is summarized by the tables below (Fields &: Olive 1999a). The table 
below represents a compilation of data using a uniform set of stellar parameters 
(effective temperature, surface gravity, and metallicity) to extract abundances. 
Results for two methods are shown (see Fields & Olive 1999a for more details). 



Table 1. Observed Pop II logarithmic slopes for Be versus Fe and O 



metal tracer 


method 


metallicity range 


Be slope 


B slope 


B/Be slope 


Fe/H 
O/H 
Fe/H 
O/H 


Balmer 
Balmcr 
IRFM 
IRFM 


-3 < \Fe/H\ < -1 
-2.5 < [O/H] < -0.5 

-3 < [Fe/H] < -1 
-2.5 < [O/H] < -0.5 


1.21 ±0.12 
1.76 ±0.28 
1.30 ±0.13 
1.38 ±0.19 


0.65 ±0.11 
1.84 ±0.58 
0.77 ±0.13 
1.35 ±0.30 


-0.18 ±0.15 
-0.81 ±0.44 
0.01 ±0.14 
0.00 ±0.17 



In standard GCRN (Meneguzzi, Audouze, & Reeves 1971), ^Be is a sec- 
ondary isotope, and is expected to have a logarithmic slope of 2 with respect to 
[Fe/H]. In the context of GCRN, B/H is also produced with a slope of 2, but 
there is an additional source for boron from neutrino spallation in supernovae 
(Hartmann 1999) which is primary (slope of 1). The [Fe/H] data shown in the 
table, however, do not reflect the expectations of standard GCRN. As a result, 
many models have been developed over the last several years to explain this 
data. For a review of these, see Ramaty (1999). 

The apparent failure of GCRN to account for the evolution of BeB vs. 
[Fe/H] rests upon the assumption that at low metallicity, [0/Fe] is constant. 
However, new data (discussed here by Garci-Lopez) from Israelian et al. (1998) 



indicates that in fact 0/Fe is not constant. If we take [0/Fe] = WQ/Fe [Fe/H] , 
then we would expect up to an additive constant (Fields & Olive 1999a) 

[Be] = 2(1 + iOo/Fc) [Fe/H] (5) 

Now for the Israelian et al. (1998) value of wq/Fc = —0.31, we would predict a 
Be slope which is consistent with the data. Galactic chemical evolution models 
with a varying 0/Fe were presented here by Fields. As we will see below, these 
models also affect the evolution of ^Li in a positive way (Fields & Olive 1999b). 

4. Two- component ^Li 

In order to test big bang nucleosynthesis, it is necessary to establish a primordial 
abundance of Li. The extraction of the primordial Li abundance is complicated 
by two factors: the GCRN production of Li and the depletion of Li in halo 
stars. If we for the moment ignore the depletion, we must still ascertain what 
fraction of the observed ^Li is primordial. In principle, we can use the abundance 
information on the other LiBeB isotopes to determine the abundance of the 
associated GCRN produced ^Li. As it turns out, the boron data is problematic 
for this purpose, as there is very likely an additional source for ^^B, namely 
i^-process nucleosynthesis in supernovae as indicated above. In Walker et al. 
(1993), and Olive & Schramm (1992), the Be data was used to set a rough 
upper bound of 20-30% to the fraction of GCRN produced '^Li on a star by star 
basis. 

There are in fact many stars for which both ^Li and ^Be have been detected. 
Using this subset of the data, one can extract the primordial abundance of Li 
in the context of a given model of GCRN. For example, a specific GCRN model, 
predicts the ratio of Li/Be as a function of [Fe/H] . Under the (plausible) assump- 
tion that all of the observed Be is GCRN produced, the Li/Be ratio would yield 
the GCRN produced ^Li and could then be subtracted from each star to give a set 
of primordial ^Li abundances. This was done in Olive & Schramm (1992) where 
it was found that the plateau was indeed lowered by approximately 0.07 dex. 
However, it should be noted that this procedure is extremely model dependent. 
The predicted Li/Be ratio in GCRN models was studied extensively in Fields, 
Olive & Schramm (1994). It was found that Li/Be can vary between 10 and 
~ 300 depending on the details of the cosmic-ray sources and propagation-e.g., 
source spectra shapes, escape pathlength magnitude and energy dependence, 
and kinematics. 

In contrast, the ^Li/^Li ratio is much better determined and far less model 
dependent since both are predominantly produced by a — a fusion rather than by 
spallation. The obvious problem however, is the paucity of ^Li data. As more ^Li 
data becomes available, it should be possible to obtain a better understanding 
of the relative contribution to ^Li from BBN and GCRN. 



5. Li Depletion 

Stellar evolution models have predicted depletion factors which differ widely, 
ranging from essentially no depletion in standard models (for stars with T ^ 
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5500 K) to a large depletion (Deliyannis et al. 1990, Charbonnel et al. 1992). 
Depletion occurs when the base of the convection zone sinks down and is exposed 
to high temperatures, ~ 2 x 10^ K for ^Li and ~ 1.65 x 10^ K for ^Li (Brown 
& Schramm 1988). In standard stellar models, the depletion of ^Li is always 
accompanied by the depletion of ^Li, though the converse is not necessarily 
true. Below, the consequences of ^Li depletion will be examined from both its 
effect on the galactic evolution of ®Li and its effect on the concordance of BBN 
and the observations of D and ^He. 

5.1. ®Li and depletion 

Model results (Fields & Olive 1999b) for ^Li vs. Fe appear in Figure 5, for an 
[0/Fe]-[Fe/H] Pop II slope wq/fc = —0.31 (as discussed above) and ujq/Fc = 
for comparison. We see that GCRN does quite well in reproducing both solar 
and Pop II ^Li when 0/Fe is allowed to evolve in Pop II. On the other hand, if 
0/Fe is constant, then the ^Li-Fe slope is steeper and the model underproduces 
the Pop II ^Li. This is consistent with what is obtained for the evolution of BeB 
(Fields & Olive 1999a,c). 
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Figure 5. The ^Li evolution as a function of [Fe/H]. Solid line: the 
"revised standard" GCRN model. Here Fe is scaled from the calculated 
O to fit the observed [0/Fe]-[Fe/H] slope. Dashed line: the GCRN 
model with Fe oc O in Pop II. The error bars on the points are 2 
sigma errors, and the spread in the points connected by lines show the 
uncertainty due to stellar parameter choices. 



It is also useful to compare the evolution of ^Li with that of Be or B. The 
problem again is the paucity of data. For HD 84937, data is available for both 
®Li and Be. By comparing the solar ratio of ^Li/Be0 = 5.9 and the value of this 
ratio for HD84937, ^Li/Be ~ 80, it appears that GCRN is consistent (though the 
observational uncertainties remain large). This data are good enough however, 
to exclude the possibility of a purely primary (linear) evolution for [Be/H] vs. 
[Fe/H] coupled with the expected linear evolution of ^Li as a primary element 
due to a — a fusion. In this case one would expect the ^Li/Be ratio to be 
constant, which is clearly not the case. 

In standard GCRN (with constant 0/Fe at low metallicities) , ^Be is a sec- 
ondary isotope, and given the linearity of pLi], one should expect that ^Li/Be 
is inversely proportional to Fe/H (i.e., to have a log slope of -1). However, if we 
take cJo/Fe 7^ 0, then we would expect the Be evolution to be governed by Eq. 
(5) (Fields & Olive 1999a,b) and 

[6Li] = (l+c^o/Fe) [Fe/H] (6) 

so that 

[6Li/Be] = -(l + u.o/Fe) [Fe/H] (7) 

The Israelian et al. (1998) value of cjo/fc = —0.31, implies a dependence which 
is consistent with the data. 

The above comparison of the data as shown in Figure 5 to the models do not 
take into account any depletion of ^Li. There is still a great deal of uncertainty 
in the amount of depletion for both ^Li and ^Li as well as the relative depletion 
factor, Dq/ Dj (Chaboyer 1994, Vauclair and Charbonnel 1995, Deliyannis et 
al. 1996, Chaboyer 1998, Pinsonneault et al. 1992, Pinsonneault et al. 1998, 
Pinsonneault 1999, Vauclair 1999). The observed lithium abundance can be 
expressed as 

Liobs = I?7(^LiBB + ^LicR) + D6(^LiBB + ^LIcr) (8) 

where the D(,j < 1 are the ^'^Li depletion factors. The two lithium components 
due to big bang and cosmic ray production, discussed in the previous section, 
are shown explicitly here. Given enough ^Li Pop II data, one could use the 
observed ^Li evolution (1) to infer '^LicR and thus ^LIbB) and (2) to measure 
^Li/Be and thereby constrain in more detail the nature of early Galactic cosmic 
rays both which would lead to a better understanding of Li depletion in general. 
In standard stellar models, Brown &: Schramm (1988) have argued that 
Dq ~ Dj with (3 ~ 60. Clearly in this case any observable depletion of Li would 
amount to the total depletion of ^Li. Hence the observation of ^Li in HD84937 
has served as a basis to limit the total amount of ^Li depletion (Steigman et 
al. 1993, Lemoine et al. 1997, Pinsonneault et al. 1998). There are however, 
many models based on diffusion and/or rotation which call for the depletion of 
^Li and Li even in hot stars. The weakest constraint comes from assuming that 
depletion occurs entirely due to mixing, so the destruction of the Li isotopes is 
the same despite the greater fragility of ^Li. Because ^Li/^Li --^ 1 in cosmic-ray 
nucleosynthesis, the observation of ^Li does exclude any model with extremely 
large ^Li depletion if one requires the preservation of the Spite plateau for Li 
up to [Fe/H] = -1.3 (Pinsonneault et al. 1998, Smith et al. 1998). However, 



barring an alternative source for the production of ^Li, the data are in fact much 
more restrictive. At the 2cj level, the model used to produce the evolutionary 
curve in Figure 5, would only allow a depletion of ^Li by 0.15 dex {Dq > 0.7); 
since Dj > Dq, this is also a lower limit to Dj. 

Further constraints on Di become available if we adopt a model which 
relates ®Li and ^Li depletion. E.g., if we use logDg = —0.19 + 1.94 log -D7 as 
discussed in Pinsonneault et al. (1998), the data in the context of the given 
model would not allow for any depletion of ^Li. However, given the observation 
al uncertainties, together with the uncertainties in the stellar parameters it is 
possible to account for some Li depletion. For example, using the Balmer line 
stellar parameters, Fields & Olive (1999a) found WQ/Pe = — 0.46±0.15. Using the 
value of -0.46, it was determined that at 2a (with respect to the ^Li data) that 
logDg > —0.32 and would still limit logDj > —0.07. Even under what most 
would assume is an extreme 0/Fe dependence of wq/Fc = —0.61, ^Li depletion is 
limited to by a factor of 3.5 and corresponds to an upper limit on the depletion 
of ^Li by 0.2 dex. This is compatible with the upper limit in Lemioine et al. 
(1997) though the argument is substantially different. 

It should be clear at this point, that improved (= more) data on ^Li in halo 
stars can have a dramatic impact on our understanding of cosmic-ray nucleosyn- 
thesis and the primordial abundance of ^Li. Coupled with improved data on the 
0/Fe ratio in these stars, we would be able to critically examine these models 
on the basis of their predictions of ^Li and ^Be. 

5.2. Constraints from BBN 

It is also in principle possible to constrain the degree of ^Li depletion from the 
concordance of the light elements produced in BBN and the observations. If 
^Li depletion were significant, then the comparison of BBN predictions to the 
observed abundances should be based on a Li abundance which is greater than 
that determined in pop II halo stars. As a result of the local minimum in the 
BBN produced '^Li abundance (at about tjiq ~ 3) the two most likely values of 
7] from ^Li (the twin peaks in the likelihood distribution of Figure 2a), are split 
farther apart as is shown in Figure 6a, where the assumed ^Li abundance is ^Li/H 
= (4.1 lb 0.1) X lO"^'^. This corresponds to a depletion factor of 0.4 dex, which 
was recently argued to be an upper limit to the ^Li depletion (Pinnsoneault et 
al. 1998, Pinnsoneault 1999). As one can see, the previous excellent agreement 
between ^Li and ^He (seen in Figure 2a) is lost. The lack of agreement is seen 
more quantitatively by comparing the two likely distributions of Figure 2b and 
6b. 

When one also considers deuterium, the concordance is further disturbed. 
For high D/H, where the agreement between D/H and "^He is good, there is 
barely any overlap with Li as can be seen in Figure 7a. The total likelihood 
function in this case, shown in Figure 7b, can be compared with that in Figure 
3b, which shows the likelihood distribution with high D/H and no Li depletion. 

Finally, one can consider the case of low D/H. Here the concordance was 
never very good. ^Li depletion causes the high-r/ peak of the Li distribution to 
move to the right (to higher values of rj), and for a depletion factor of 0.4 dex 
as is assumed in Figure 8, the peak of the ^Li distribution lies at values of rj 
larger than that predicted by low D/H. Clearly for a smaller ^Li depletion factor, 
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good agreement between ^Li and D/H can be achieved. In this case, ^He is still 
somewhat problematic. 




Figure 6. (a) - As in Figure 2a, where the assumed primordial ^Li 
abundance has been increased to take into account a possible 0.4 dex 
of depletion; (b) - As in Figure 2b, with the increased primordial ^Li. 




Figure 7. (a) - As in Figure 6a, with the addition of the likelihood 
distribution for D/H assuming high D/H; (b) - The total likelihood 
distribution with high D/H. 



6. Summary 

It is apparent that we have a general understanding of the production mecha- 
nisms of the LiBeB elements. Their (relative) inferior position on the abundance 
chart indicates that they are not produced in the normal course of stellar evolu- 
tion and as is well known the mass gaps at A = 5 and 8 prevent them from being 
produced in sufficient abundance in the big bang. Cosmic-ray nucleosynthesis 
has been shown to work reasonably well in accounting for the LiBeB abundances. 
However, several questions remain. Among them is the challenge to separate the 
BBN and GCRN components of ^Li. This is crucial for testing the concordance 
of BBN theory with observations. The depletion of ^Li is another complication 
that must be resolved. It was argued above, that ^Li may hold the key to re- 
solving both of these problems. The similar processes (namely a — a fusion) 
which produce both Li isotopes could allow for a direct determination of GCRN 
produced ''^Li. Due to its fragility, ^Li, if observed in more halo stars could also 
play a key role in pinning down the Li depletion factor. In addition, the com- 
parison with BeB, would be extremely useful in distinguishing between primary 
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Figure 8. (a) - As in Figure 6a, with the addition of the hkehhood 
distribution for D/H assuming low D/H; (b) - The total likelihood 
distribution with low D/H. The dashed curve is the distribution from 
Figure 7b. 



and secondary models of cosmic-ray nucleosynthesis. While great advances have 
occurred in obtaining LiBeB data, open questions can only be resolved with new 
data on these elements. 
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